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ABSTRACT 


A  high  resolution,  20-level,  primitive  equation  (PE)  model  of  the  California  coastal 
region  is  initialized  using  temperature  data  acquired  during  the  first  Naval  Postgraduate 
School  California  Transition  Zone  (CIZ)  cruise  from  6  to  12  July  1988  for  the  purpose 
of  diagnosing  the  three-dimensional  dynamically  balanced  flow  field  for  the  region.  The 
major  feature  in  the  region  during  the  cruise  period  was  a  strong  meandering  jet  which 
flowed  equatorward  and  offshore,  oriented  in  a  northeast  to  southwest  direction.  The 
quality  of  data  acquired  during  cruise  CTZRI  is  sufficiently  high  to  enable  diagnosis  of 
the  horizontal  (baroclinic)  and  vertical  velocity  field  using  the  numerical  model.  The 
measured  (ADCP)  currents  showed  more  details  of  the  flow  at  deeper  depths  than  the 
model,  which  showed  a  broader  baroclinic  flow  at  depth  and  a  level  of  no  motion  near 
300  m.  The  maximum  surface  velocities  from  model  and  ADCP  cross-sections  agreed 
to  within  10  cm  sec-1  except  for  leg  G,  in  which  ADCP  velocity  was  greater  than  the 
model  velocity  by  about  30  cm  sec-1.  The  sign  of  vertical  velocity  agrees  very  well  with 
independent  estimates  made  from  bio-optical  data,  however,  the  magnitude  calculated 
by  the  model  is  30  to  60  times  larger  th  .n  that  estimated  from  the  bio-optical  observa¬ 
tions. 
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I.  INTRODUCTION 


A.  THE  CALIFORNIA  CURRENT  SYSTEM 

The  California  Current  is  the  equatonvard  leg  of  the  anticyclonic  gyre  covering  the 
majority  of  the  North  Pacific  between  5’  to  10“  N  and  45*  to  50’  N  (Tchernia,  I9S0).  The 
core  of  the  mean  California  Current  is  located  approximately  100  to  200  km  offshore  in 
the  upper  200  m.  Current  velocities  range  from  5  cm  sec-1  to  14  cm  sec*1  depending  on 
the  season  and  location  (Chclton,  1984).  Chelton  calculated  these  velocities  using 
geostrophv,  with  500  m  as  the  level  of  no  motion.  Over  30  years  of  hydrographic  survey 
data  from  California  Cooperative  Oceanic  Fisheries  Investigations  between  San 
Francisco  and  Point  Conception  were  used  in  his  study.  In  addition  to  the  equatorward 
flowing  California  Current  there  exists  the  poleward  flowing,  subsurface  California 
Undercurrent;  the  poleward  flowing,  surface  Davidson  Current,  found  north  of  Point 
Conception  during  fall  and  winter;  and  the  poleward  flowing,  surface  Southern 
California  Countercurrent,  found  south  of  Point  Conception  in  the  California  Bight. 
These  four  currents  make  up  the  California  Current  System  (CCS).  Hickey  (1979)  pro¬ 
vides  a  thorough  description  of  the  seasonal  variation  of  the  currents  in  the  CCS. 

The  CCS  is  important  to  the  climate  and  economy  of  the  west  coast  states.  Up- 
welling  events  bring  up  cold,  nutrient-rich  subsurface  water  to  the  ocean  surface,  pro¬ 
viding  nourishment  to  the  phytoplankton  in  the  euphotic  zone  (Brink,  1983;  Buyer, 
1983;  Parsons  ei  a  ,  1977).  The  cool  sea-surface  temperatures  adjacent  to  the  coast 
engenders  mild  weather  conditions  and  the  ocean  provides  moisture  for  the  surface 
boundary  layer,  promoting  cloud  formation  and  fog  (Foster,  19S9).  Additionally,  sur¬ 
face  salinity  fronts  bounding  the  offshore  edge  of  a  meandering  jet;  temperature  fronts; 
and  sharp,  near-surface  density  fronts  associated  with  upwelling  affect  sonar  perform¬ 
ance  in  fleet  operating  areas.  The  continuing  study  of  the  CCS,  representing  a  typical 
eastern  boundary  current,  is  necessary  for  the  promotion  of  both  military  and  civilian 
interests. 

Interesting  features  of  the  CCS  are  the  cold  filaments  occurring  off  the  west  coast 
of  the  United  States,  as  far  north  as  Vancouver  Island,  British  Columbia  and  as  far  south 
as  Point  Conception,  California  (Ikeda  and  Emery,  1984).  These  filaments  (also  known 
as  tongues,  jets,  or  squirts)  consist  of  cold,  salty  inshore  water  which  has  been  carried 
offshore,  beyond  the  continental  shelf  (Brink,  1983;  CTZ  Group,  1988;  Ikeda  and  Emery, 
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1984).  Summer  infrared  satellite  images  from  the  West  Coasts  of  California  and  Oregon 
have  shown  these  filaments  extending  westward  from  the  coast  (Ikeda  and  Emery,  1984; 
CTZ  Group,  1988;  Ramp  and  Jessen,  1988;  Stanton  and  Stockel,  1988). 

The  purpose  of  the  Coastal  Transition  Zone  (CTZ)  program  is  to  investigate  the 
cold  filaments  associated  with  the  CCS.  Four  major  goals  of  the  CTZ  program  are  to 
understand  the  processes  involved  in  the  interaction  of  the  filaments  with  open  ocean 
water;  determine  the  effect  of  the  filaments  on  the  biology  and  meteorology  of  the 
coastal  zone;  create  a  series  of  maps  of  the  current,  hydrographic,  nutrient,  and  particle 
fields;  and  discover,  through  use  of  numerical  models,  what  driving  forces,  dynamics, 
topography,  and  stratification  is  necessary  to  duplicate  the  observed  characteristics  of 
these  filaments  (CTZ  Group,  1988). 

B.  OBJECTIVE 

The  objective  of  this  thesis  is  to  diagnose  the  velocity  fields  from  acquired  temper¬ 
ature  data  in  the  region  off  Point  Arena  from  approximately  37‘10'  to  39*2 l'N  and 
123*50'  to  1 26*40' W,  which  covers  roughly  one  meander  of  the  California  Current.  Areas 
of  convergence  and  divergence  and  their  relation  with  the  diagnosed  vertical  velocity 
fields  will  be  examined.  Knowledge  of  the  velocity  fields  is  important  to  tracer  studies, 
and  the  vertical  velocity  field  is  of  continuing  interest  to  biologists  and  the  local  fisheries 
industry. 
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II.  THE  MODEL 


The  numerical  model  is  a  20-level,  primitive  equation  model  of  a  baroclinic  ocean 
on  an  f-plane,  with  hydrostatic,  Boussinesq,  and  rigid  lid  approximations.  The  space- 
staggered  B-scheme  (Messinger  and  Arakawa,  1982)  is  used  for  horizontal  finite  difler- 
encing,  and  a  sigma  coordinate  system  defines  the  vertical  dimension.  The  boundary 
conditions  at  the  sea  surface  consist  of  no  wind  stress  and  no  heat  flux.  The  lateral 
boundaries  of  the  model  domain  are  treated  as  open  boundaries  using  the  method  of 
Ross  and  Orlanski  (19S3).  The  total  depth  of  the  ocean  is  4000  m,  with  model  levels 
listed  in  Table  1.  The  values  of  constants  used  in  the  model  are  listed  in  Table  2.  No 
attempt  is  made  to  model  the  depth  averaged  part  of  the  currents  which  are  taken  to  be 
zero. 


Table  1,  K-LEVELS  AND  CORRESPONDING  DEPTHS  IN  ME¬ 


TERS  USED  IN  MODEL 


k  level 

z(k)  (meters) 

k  level 

z(k)  (meters) 

1 

5 

11 

1105 

a 

4. 

25 

12 

1325 

3 

65 

13 

1565 

4 

125 

14 

1825 

5 

205 

15 

2105 

6 

305 

16 

2405 

7 

425 

17 

2725 

8 

565 

18 

3065 

9 

725 

19 

3425 

10 

905 

20 

3805 
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Table  2.  CONSTANTS  USED  IN  MODEL 


Constant 

Value 

Name 

Q. 

2n  day1 

earth  rotation  rate 

DTAU 

600  s 

timestep 

dx 

4  x  105  cm 

meridional  grid  spacing 

dv 

4  x  10J  cm 

zonal  grid  spacing 

D 

4  x  105  cm 

total  ocean  depth 

LAT 

38.25'  N 

reference  latitude 

f 

9.00  x  10*s  sec-1 

Coriolis  parameter 

g 

980  cm  sec-2 

acceleration  of  gravity 

mm 

BB 

biharmonic  momentum  diffusion 
coefficient 

mm 

4  x  10p  cm4  sec-1 

biharmonic  heat  diffusion  coeffi¬ 
cient 

KM 

0.5  cm 2  sec*1 

vertical  eddy  viscosity 

K„ 

0.5  cm2  sec-1 

vertical  eddy  conductivity 

a 

2.877  x  10*4  C*1 

thermal  expansion  coefficient 

The  chart  of  the  1988  CTZ  station  grid  is  shown  in  Figure  1.  The  station  grid  ex¬ 
tends  approximately  120  nm  offshore  and  110  nm  alongshore,  roughly  parallel  to  the 
California  coast  and  southwest  of  Point  Arena.  The  alongshore  sections  are  approxi¬ 
mately  20  nm  apart  and  have  individual  stations  approximately  13  to  14  nm  apart.  The 
model  domain  covers  the  area  from  37T0'  to  39’21'N  and  123*50'  to  1 26*40' W.  The 
chart  of  the  model  domain  with  included  stations  is  shown  in  Figure  2.  The  model  do¬ 
main  excludes  stations  Al,  A2,  A3,  A4,  A5,  B8,  and  B9.  The  actual  topography  in  the 
model  domain  is  shown  in  Figure  3.  The  effect  of  variable  tomography  and  depth  aver¬ 
aged  currents,  as  well  as  wind  forcing  and  surface  heat  fluxes  is  considered  a  separate 
study. 
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Figure  1,  CTZ  Station  Grid  for  1988:  The  irregular  alongshore  leg  paralleling  the 
coast  is  leg  A,  followed  by  legs  B,  C,  D,  E,  F,  and  G.  Stations  are  num¬ 
bered  from  the  top  to  the  bottom  of  each  leg.  Leg  A  has  12  stations,  legs 
B  and  C  have  nine  stations  each,  legs  D,  E,  and  F  have  ten  stations  each, 
and  leg  G  has  five  stations  {6  to  10)  (from  Huyer  et  al,  1990). 


CHART  OF  MODEL  DOMAIN  AND  INCLUDED  STATIONS 


Figure  2.  Chart  of  Model  Domain  with  Included  Stations:  Starting  from  left  tc 
right,  the  alongshore  sections  are  G,  F,  E,  D,  C,  B,  and  A.  Stations  A1 
A2,  A3,  A4,  A5,  B8  and  B9  are  excluded  from  the  model  domain. 
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Figure  3.  Actual  Topography  in  the  Model  Domain:  The  box  shows  the  area  of 
the  ocean  off  Point  Arena  in  the  model  domain.  Depths  are  in  meters. 


III.  DATA  ACQUISITION  AND  PROCESSING 


A.  FIELD  ACQUISITION 

The  three  Naval  Postgraduate  School  cruises  were  CTZR1  (6  to  12  July  1988), 
CTZR2  (13  to  IS  July  1988),  and  CTZ2  (22  to  26  July  1988),  all  aboard  R/V  POINT 
SUR.  Both  acoustic  doppler  current  profiler  (ADCP)  data  and  conductivity,  temper¬ 
ature  and  depth  (CTD)  data  were  acquired.  This  study  will  concentrate  predominantly 
on  the  results  from  cruise  CTZR1. 

CTD  measurements  onboard  R/V  POINT  SUR  were  made  with  an  NBIS  CTD 
system.  Temperature  measurement  accuracy  was  better  than  ±0.01  C  and  salinity 
measurement  accuracy  was  about  ±  0.003  psu.  An  ADCP  made  by  RD  Instruments 
was  operated  continuously  onboard  R/V  POINT  SUR.  The  transducer  operated  at  150 
kHz  with  a  4  meter  pulse  length  and  a  nominal  range  of  350  meters.  Returning  echoes 
were  sampled  selectively  in  4  meter  bins  to  produce  a  vertically-averaged,  ship-relative 
velocity.  Relative  velocity  profiles  were  obtained  once  every  three  minutes.  Twenty-five 
meters  was  the  shallowest  level  to  yield  reliable  data  (Huyer  et  al,  1990). 

Figure  4  is  an  infrared  Advanced  Very  High  Resolution  Radiometer  (AVHRR)  im¬ 
age  from  the  National  Oceanic  and  Atmospheric  Administration  (NOAA)  NOAA-9 
satellite,  taken  09  July  1988  of  the  ocean  off  the  coast  of  northern  California.  Cape 
Mendocino,  Point  Arena,  and  Point  Reyes  can  be  seen.  The  large  filament  oriented 
northeast  to  southwest  off  of  Point  Arena  is  the  filament  under  study. 

An  infrared  image  of  the  same  area  seven  days  later,  taken  during  the  middle  of 
cruise  CTZR2,  is  shown  in  Figure  5. 

Two  color-enhanced  images  of  the  ocean  off  Point  Arena  are  shown  in  Figure  6  and 
Figure  7.  These  two  images  were  taken  before  and  after  cruise  CTZ2. 

The  four  images  show  the  spatial  and  temporal  evolution  of  the  complex  sea  surface 
temperature  field.  The  prevailing  northwesterly  winds  in  the  summer  cause  Ekman 
transport  of  coastal  waters  offshore,  which  induces  the  upwelling  of  cooler  subsurface 
water.  This  cool  subsurface  upwelled  water  is  indicated  by  the  lighter  shades  in  the  black 
and  white  images  and  the  blues  and  purples  in  the  color-enhanced  images.  Two  large 
meanders  of  the  California  Current  can  be  seen  in  the  09  July,  16  July,  and  19  July  im¬ 
ages,  with  the  filament  under  study  separating  the  two  meanders.  In  the  black  and  white 
images,  eddy-like  features  can  be  seen  offshore  south  of  Cape  Mendocino,  Point  Arena, 
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Figure  4.  AVHRR  Infrared  Image  of  the  Ocean  Off  the  Coast  of  Northern 
California,  09  July  1988:  Darker  shades  denote  warmer  temperatures. 
Lighter  shades  denote  cooler  temperatures.  The  latitude  coordinates 
start  at  37‘N  at  the  bottom  left  corner  and  increase  to  41°X.  The  longi¬ 
tude  coordinates  start  at  123‘W  at  the  bottom  right  corner  and  increase 
to  12S'W.  Satellite  image  provided  courtesy  of  Prof.  Steven  R.  Ramp 
and  Mr.  Paul  Jesscn,  Oceanography  Department,  Naval  Postgraduate 
School. 
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Figure  5.  AVHRR  Infrared  Image  of  the  Ocean  Off  the  Coast  of  Northern 
California,  16  July  1988:  Darker  shades  denote  warmer  temperatures. 
Lighter  shades  denote  cooler  temperatures.  The  latitude  coordinates 
start  at  37'N  at  the  bottom  left  corner  and  increase  to  41CN.  The  longi¬ 


tude  coordinates  start  at  123:\Y  at  the  bottom  right  corner  and  increase 


to  128  \V.  Satellite  image  provided  courtesy  of  Prof.  Steven  R.  Ramp 
and  Mr.  Paul  Jcsscn,  Oceanography  Department,  Naval  Postgraduate 


School. 
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Figure  6.  Color- Enhanced  AVHRR  Infrared  Image  of  the  Ocean  Off  the  Coast  of 
Northern  California,  19  July  1988:  Purples  and  blues  denote  cooler 

temperatures.  Violets  and  reds  denote  warmer  temperatures.  Satellite 
image  provided  courtesy  of  Prof.  Timothy  P.  Stanton  and  Mr.  James 
Stockel,  Oceanography  Department,  Naval  Postgraduate  School. 

and  Point  Reyes.  These  satellite  images  help  to  give  an  appreciation  of  the  complex  way 
in  which  the  California  Current  interacts  with  the  coastal  waters  to  produce  the  fila¬ 
ments  and  eddies  evident  in  these  images. 
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Figure  7.  Color-Enhanced  AVHRR  Infrared  Image  of  the  Ocean  Off  the  Coast  of 
Northern  California.  27  July  1988:  Purples  and  blues  denote  cooler 

temperatures.  Violets  and  reds  denote  warmer  temperatures.  Satellite 
image  provided  courtesy  of  Prof.  Timothy  P.  Stanton  and  Mr.  James 
Stockcl,  Oceanography  Department,  Naval  Postgraduate  School. 

B.  DATA  PROCESSING 

The  acquired  temperature  data  was  objectively  analyzed  using  optimum  interpo¬ 
lation  as  first  detailed  by  Gandin  (1965)  for  meteorological  use  and  later  applied  to 
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oceanographic  data  by  Bretherton  et  al„  (i976)  and  Carter  and  Robinson  (1987).  The 
method  has  the  following  features  (Stanton  and  Stockel,  1988): 

Provides  a  consistent  method  to  represent  sparsely  sampled  fields; 

Forms  an  optimal  estimate  of  an  n-dimensional  variable  at  gridded  intervals  based 
on  a  measured  and/or  modeled  correlation  function; 

Provides  a  measure  of  the  error  of  the  estimated  variable  at  each  grid  point; 

Allows  straightforward  interfacing  to  models  and  calculations  involving  field  differ¬ 
ences  and  integration,  and  'standardizes'  criteria  for  graphical  representation  of 
mapped  data  fields. 

Processing  the  acquired  data  for  initialization  of  the  20-level  model,  which  carries 
only  temperature  as  a  scalar,  proceeded  as  follows.  The  mean  temperature  and  mean 
salinity  profiles  for  each  of  the  three  cruises  was  calculated.  The  acquired  temperature 
data,  T(z),  was  used  with  the  domain  mean  salinity  profile,  S(z),  in  an  iterative  scheme 
to  produce  an  apparent  temperature  profile,  T.(z),  for  input  into  the  model.  The  full 
equation  of  state  was  used  in  the  iterative  scheme  to  derive  Ta{z).  The  density  calculated 
from  the  apparent  temperature  and  domain  mean  salinity,  p{Tc(z),  S(z)),  was  the  same 
as  the  in  situ  density,  p(T(z),  S(z)),  measured  at  each  station.  The  apparent  temper¬ 
atures  were  vertically  averaged  in  four-meter  bins  at  5  m,  25  m,  65  m,  125  m,  205  m,  305 
m,  and  425  m  (the  first  seven  k-levels  in  the  model)  in  the  optimum  interpolation  pro¬ 
gram  for  cruises  CTZR1  and  CTZ2.  The  apparent  temperatures  were  vertically  averaged 
in  four-meter  bins  at  5  m,  25  m,  65  m,  125  m,  and  205  m  for  cruise  CTZR2. 

The  spatial  m  an  temperature  profile  was  extended  to  4000  m  by  fitting  an  expo¬ 
nential  tail  to  the  last  two  points  of  the  calculated  mean  temperature  profile.  This  ex¬ 
tended  mean  temperature  profile  was  used  to  calculate  the  first  three  dynamical  vertical 
modes  T2{k )  ,  and  T^k)  for  temperature.  The  objective  was  to  extrapohtc  the 
perturbation  temperature  down  to  4000  m  using  modal  profiles  with  weights  determined 
according  to  how  well  the  data  in  the  upper  425  m  (205  m  for  CTZR2)  fit  the  respective 
mode.  The  procedure,  as  outlined  below,  is  described  in  detail  in  Rienecker  et  ah, 
(1987). 

The  initial  assumption  made  was  that 

TVJ,k)=T(i,j,k)-T(k),  {1} 

i,  j,  k,  are  the  gridpoint  indices  for  the  x,  y,  and  z  directions,  respectively, 

T '(/ ,  j ,  k  )is  the  perturbation  part  of  the  temperature  at  the  gridpoint, 
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T(i  J  ,k)  is  the  objectively  analyzed  temperature  at  the  gridpoint,  and 
T(k)  is  the  mean  temperature  at  level  k. 

The  perturbation  temperature  T '(/ ,  j ,  k )  is  assumed  to  be  a  linear  combination  of 
the  first  three  dynamical  vertical  modes,  such  that 

T '(/  J  ,k )  =  «,(/  J)T}(k)  +  a2(i  ,j)T2(k)  +  a3(i  ,j)T3(k),  {2} 

ai  ~  ai(i  J)  is  the  weighting  given  to  7j(A) , 
a2  =  a2{i  J)  is  the  weighting  given  to  T2(k) ,  and 
ai  -  ai(‘  J)  is  the  weighting  given  to  T2{k) . 

Solutions  for  ax,  a2,  and  a}  are  desired. 

Four  equations  are  derived  from  equation  {2}: 

<  T '(/  ,j  ,k)>  =a1<Tl(k)>  +a2<T2(k)>  +a3<T3(k)>,  {3} 

<T,(k)T'(i,j,k)>  =a,<  T1(k)Tl(k)>  +  a2<  Tl(k)T2(k)>  +  a3<  Tl(k)T3(k)> ,  {4} 
<T2(k)T'(iJ,k)>  =  <  T2(k)T]{k)>  +  a2<T2(k)T2(k)>  +  a3<T2(k)T3(k)> ,  {5} 

<T3(k)TVJ,k)>  =al<r3(A)r1(A)>  +a2<T3(k)T2(k)>  +  a3<T3{k)T3(k)> ,  {0} 

i  ■> 

and  define  the  bracketed  terms,  <Q>  s  ^gA2(A)/£Az(A),  where  A z{k)  is  the  vertical 

i  i 

difference  between  adjacent  k  levels  in  the  model,  Thus,  <  Q  >  is  the  vertical  average 
of  Q  over  the  top  seven  levels  of  the  model  which  contain  all  the  analyzed  observations. 

The  system  of  four  equations  in  three  unknowns  (a„  a2,  a3)  can  be  solved  by  linear 
least  squares  methods.  After  a„  a2 ,  and  a3  are  found  for  each  point  in  the  domain  of  the 
model,  equation  {2}  can  then  be  used  to  find  T'(i,j,k)  for  k  =  8  to  k  =  20.  For 
CTZR2,  the  procedure  is  similar,  except  only  five  levels  of  data  are  available,  and  the 
extrapolation  is  used  to  define  T'(i  ,j ,  k )  for  k  =  6  to  k  =  20. 

Solutions  to  equations  {3}  through  {6}  above  yielded  values  of  a,  which  produced 
negative  temperatures  for  T(i  ,j  ,k )  at  k  levels  9  through  14  (725  m  to  1825  m),  which 
caused  the  mod  ;i  to  become  unstable  at  42  time  steps  (0.29  day).  Apparently,  the  CTD 
casts  to  only  500  meters  were  not  sufficient  to  adequately  define  the  first  dynamical 
mode  to  4000  m.  Accordingly,  the  new  system  of  equations  used  to  solve  for  a2  and  a3 
were  the  following: 

<  T '(/  ,j  ,k)  >  =a2<T2(k)>  +a3<T3(k)> ,  {7} 
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<Uk)T’uj,k)>  -fl2<r,(*)r#)>  +a3<uk)T3(k)>, 


{8} 


<T2(k)T'(iJ,k)>  =  a2<T2{',)T2{k)>  +  a3<T2{k)T3(k)> ,  {9} 

<T3{k)T'(i,j,k)>  =  <  T3(k)T2(k) >+a3<  T3(k)T3(k) > .  {10} 

The  resulting  solutions  for  a2  and  Aj  provided  positive  extrapolated  temperatures  in 
the  model  domain;  however,  horizontal  temperature  fluctuations  on  scales  less  than  the 
observation  scale  (about  25  km)  were  apparently  hitroduced  during  the  objective  analy¬ 
sis  scheme.  This  was  discovered  after  contouring  the  extrapolated  input  temperatures 
as  seen  in  Figure  8. 

If  these  small  scales  are  in  fact  physical,  then  their  vertical  scales  should  also  be 
small  (much  smaller  than  that  of  dynamical  modes  2  and  3).  Since  the  extrapolation 
scheme  uses  modes  2  and  3  to  extrapolate  the  perturbation  temperatures  down  to  4000 
m,  the  assumption  is  that  the  horizontal  scales  are  those  of  the  deep  modes,  which  are 
larger  than  25  km.  Thus  the  T '(/ ,  j ,  k )  fields  were  smoothed  to  eliminate  the  small 
scales.  The  following  smoother  was  used: 

TS~-~{4T0+T+  +  TX),  {11} 

T0  =  T'(i  ,j ,  k ),  the  perturbation  temperature  at  each  domain  gridpoint, 

7  -  -  TV  +  1  J,k)+  TV  -  1  J ,  k)  +  TV  J+Uk)+  TV  J~  1 .  *), 
T'=T'(i+lJ+\,k)  +  T'(i-l,j+l,k)  +  TV-lJ-l,k)  +  T'(i+l,j-l,k), 

and  smoothing  was  accomplished  by  setting  T'(i ,j,k)s  T“  (double  smoothing).  An 
example  of  the  resulting  smoothed  field  is  shown  in  Figure  9. 

An  example  of  the  extrapolation  of  the  apparent  temperature  to  deep  depths  is  now’ 
given.  Two  points  from  the  CTZR1  five  meter  plot  of  apparent  temperature,  (52,  16) 
and  (24S,  220),  will  be  used  as  examples. 

Point  (52,  16)  is  in  a  warm  region  located  in  the  southwestern  corner  of  the  grid. 
Point  (248,  220)  is  in  a  cold  eddy  feature  in  the  northeastern  corner  of  the  grid,  off  Point 
Arena.  The  locations  of  the  two  points  are  shown  m  Figure  10. 

The  modal  temperature  values  for  CTZR1  are  shown  in  Figure  11.  The  third  mode 
has  the  highest  amplitude  (about  44‘C)  followed  by  the  second  mode  (maximum  ampli¬ 
tude  about  18”C)  and  the  first  mode  (maximum  amplitude  about  3'C).  The  first  zero 
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Figure  8.  Contour  Plot  of  Extrapolated  Apparent  Temperature  at  5o5  m  (un¬ 
smoothed):  The  numerals  on  the  abscissa  and  ordinate  refer  to  dis¬ 
tances  in  kilometers.  The  model  domain  only  includes  the  i  indices  from 
0  to  256  and  the  j  indices  from  1  to  240.  The  irregular  contours  in  the 
northeast  corner  are  the  zero  contours  which  designate  the  coastline 
(Point  Arena  is  outlined). 

crossing  of  the  third  mode  is  at  about  425  m  and  the  second  zero  crossing  is  at  about 
1250  m.  The  zero  crossing  of  the  second  mode  is  at  about  800  m. 
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Figure  9.  Contour  of  Extrapolated  Apparent  Temperature  at  565  m 
(smoothed):  The  numerals  on  the  abscissa  and  ordinate  refer  to  dis¬ 


tances  in  kilometers.  The  model  domain  only  includes  the  i  indices  from 


0  to  256  and  the  j  indices  from  0  to  240.  The  irregular  contours  in  the 
northeast  corner  are  the  zero  contours  which  designate  the  coastline 


(Point  Arena  is  outlineu). 


The  plots  of  the  a2  and  a%  values  for  cruise  CTZR1  are  shown  in  Figure  12  and 
Figure  13,  respectively.  The  fields  are  clearly  related  to  the  synoptic  scale  temperature 
patterns  shown  in  Figure  8  and  Figure  9. 
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Figure  10.  Plot  of  Apparent  Temperature  at  Five  Meters  Depth,  Cruise  CTZR1, 
06-12  July  1988:  The  heavily  outlined  area  in  the  northeast  comer  is 
the  coastline.  The  two  points  of  interest  are  marked  with  bullseyes. 
The  model  domain  only  includes  the  area  from  x  =  0  to  x  =  256  km 
and  y  =  0  to  y  =  24C  im, 


Note  that  for  point  (52, 16),  the  a3  value  is  positive  and  the  a2  value  is  negative,  with 
|a3|  >  \a2\,  This  implies  that  the  perturbation  temperature  at  (52,  16)  is  strongly  in¬ 
fluenced  by  the  third  mode  and  will  be  positive  until  the  first  zero  crossing  of  mode  three. 
After  the  first  zero  crossing,  the  perturbation  temperature  values  will  be  negative  until 
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the  second  zero  crossing.  After  the  second  zero  crossing,  the  perturbation  temperatures 
will  be  positive  to  4000  m. 
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50.0 


PLOT  OP  R2  VfiLUES 
CRUISE:  CTZR1,  6-12  JUL  88 


Figure  12.  Plot  of  A2  Values,  Cruise  CTZR1,  06-12  July  1988 

For  point  (248,  220),  the  az  value  is  negative  and  the  a2  value  is  positive,  with 
I  a3 1  >  |  a2 1 .  This  implies  that  the  perturbation  temperatures  at  (248,  220)  will  be  nega¬ 
tive  until  the  first  zero  crossing  of  the  third  mode,  positive  until  the  second  zero  crossing, 
and  negative  to  4000  m  after  the  second  zero  crossing. 

The  plot  of  perturbation  apparent  temperature  versus  depth  for  point  (52,  16)  is 
given  in  Figure  14.  The  plot  of  perturbation  apparent  temperature  versus  depth  for 
point  (248,  220)  is  given  in  Figure  15.  Total  and  mean  apparent  temperature  profiles 
for  points  (52,  16)  and  (248,  220)  are  plotted  in  Figure  16  and  Figure  17. 

The  point  (248,  220)  is  an  example  of  a  location  where  the  objective  analysis  resulted 
in  a  hydrostatically  unstable  profile  of  apparent  temperature.  Any  hydrostatic  instability 
was  dealt  with  by  subjecting  the  respective  point  to  a  convective  adjustment  which  cor- 
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PLOT  OF  R3  VRLUES 
CRUISE  CTZR1 ,  6-12  JUL  88 


Figure  13.  Plot  of  A3  Values,  Cruise  CTZRl,  06-12  July  1988 


rected  the  instability.  The  result  of  this  convective  adjustment  for  point  (248,  220)  is 
shown  in  Figure  18.  The  temperature  profile  is  now  hydrostatically  stable.  Together, 
these  figures  show  that  the  vertical  extrapolation  procedure  was  effective  in  extending 
the  perturbation  temperatures  to  4000  m. 

The  need  for  a  convective  adjustment  for  some  of  the  data  points  is  attributable  to 
the  correlation  function  used  in  the  objective  analysis  scheme.  The  correlation  function 
is  an  even,  smooth,  monotonic  function  with  a  zero  crossing  in  the  horizontal  space  di¬ 
mension  modeled  from  observed  data.  An  instability  in  an  objectively  analyzed  profile 
can  occur  if  the  conelation  function  encounters  a  quickly  changing  monotonic  profile 
in  the  negative  portion  of  the  correlation  function.  The  large  negative  contribution 
created  by  the  product  of  the  correlation  function  and  the  value  of  the  profile  at  a  spe- 
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Figure  14.  Perturbation  Apparent  Temperature  Versus  Depth  for  Point  (52,  16) 

cific  depth  may  be  enough  to  cause  the  once  stable  profile  to  develop  spurious 
hydrostatic  instability.  For  cruise  CTZR1,  a  plot  of  the  points  requiring  a  comective 
adjustment  is  given  in  Figure  19.  Most  of  the  convective  adjustment  occurred  north  and 
south  of  the  jet.  Cruise  CTZR2  required  convective  adjustment  at  all  data  points. 
Cruise  CTZ2  required  convective  adjustment  on  either  side  of  the  jet  ( the  jet  had  reori¬ 
ented  itself  from  a  northeast-southwest  direction  to  a  nearly  north-south  direction)  as 


PER  TURBOT  ION  FIPPRRCHT  TEMPERATURE  (C) 
VERSUS  DEPTH  (H) 

CRUISE  CTZR1 
06-12  JULY  1988 


Figure  15.  Perturbation  Apparent  Temperature  versus  Depth  for  Point  (248, 220) 

seen  in  Figure  20.  The  convective  adjustment  at  point  (248,  220)  for  both  CTZR2  and 
CTZ2  are  shown  in  Figure  21  and  Figure  22. 
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Figure  16.  Plot  of  Total  and  Mean  Apparent  Temperature  for  Point  (52, 
16):  Total  apparent  temperature  is  plotted  using  squares  and  a  solid 
line.  Mean  apparent  temperature  is  plotted  using  circles  and  a  dashed 
line.  The  squares  and  circles  are  plotted  at  the  k-level  depths  of  the 
model. 
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Figure  17,  Plot  of  Total  and  Mean  Apparent  Temperature  for  Point  (248, 
220):  Total  apparent  temperature  is  plotted  using  squares  and  a  solid 
line.  Mean  apparent  temperature  is  plotted  using  circles  and  a  dashed 
line.  The  squares  and  circles  are  plotted  at  the  k-level  depths  of  the 


model. 
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Depth  (m) 


Figure  19.  Areas  Requiring  Convective  Adjustment  for  CTZRl:  The  hatched 
areas  denote  the  points  which  required  convective  adjustment  before 
initiating  the  model  run. 


Distance  (km) 


tadj  at  day  0  depth  0  ci  =  1.00 
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Figure  21.  Plot  of  Total  Apparent  Temperature  for  Point  (248,  220)  after 
Convective  Adjustment,  Cruise  CTZR2:  The  solid  line  and  squares 
denote  the  hydrostatically  unstable  profile  after  extrapolation  and  be¬ 
fore  convective  adjustment.  The  dashed  line  and  circles  denote  the 
hydrostatically  stable  profile  after  convective  adjustment. 
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depth  M 


Total  Apparent  Temperature  (C) 
Versus  Depth  (m) 

Cruise  CTZ2 
22—26  July  1988 


Figure  22.  Plot  of  Total  Apparent  Temperature  for  Point  (248,  220)  after 
Convective  Adjustment,  Cruise  CTZ2:  The  solid  line  and  squares  de¬ 
note  the  hydrostatically  unstable  profile  after  extrapolation  and  before 
convective  adjustment.  The  dashed  line  and  circles  denote  the 
hydrostatically  stable  profile  after  convective  adjustment. 
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IV.  MODEL  RESULTS 


The  primitive  equation  model  was  initialized  with  the  apparent  temperature  fields 
described  above  and  geostrophically  balanced  currents  having  a  zero  vertical  average. 
The  model  was  then  integrated  forward  in  time  for  30  hours  ( t  P /'■)  with  the  temper¬ 
ature  field  held  fixed.  The  currents  underwent  adjustments  to  a  new  state  of  balance  in 
which  nonlinearity  and  eddy  diffusion  effects  are  not  entirely  negligible  (  as  they  are  in 
exact  geostrophic  balance).  The  resulting  density  diagnostic  currents,  which  have  an 
important  vertical  component,  are  the  main  subject  of  analysis  and  verification. 

A.  PRESSURE,  TEMPERATURE,  AND  VELOCITY  FIELDS 

The  discussion  of  model  results  will  be  limited  to  the  model  output  from  CTZR1 
data  and  the  vertical  velocity  output  at  time  30  hours  for  both  CTZR2  and  CTZ2  only. 
Fuller  treatment  of  cruises  CTZR2  and  CTZ2  are  considered  to  be  separate  studies. 

The  apparent  temperature  fields  at  50  m,  100  m,  and  200  m  at  time  zero  are  shown 
in  Figure  23,  Figure  24,  and  Figure  25.  Ver  ;cal  continuity  is  evident  in  the  apparent 
temperature  fields.  In  each  level  a  strong  thermal  gradient  oriented  southeast  to  north¬ 
west  appears,  with  a  strong  jet  dividing  the  domain  into  two  distinct  parts.  This  strong 
apparent  temperature  gradient  (VT «  l*/8  km  a;  50  m  and  1*/16  km  at  200  m)  is  the 
driving  force  for  the  strong  initial  currents  in  the  upper  ocean.  These  currents  are  com¬ 
puted  from  the  thermal  wind  relation  and  with  the  assumption  that  the  vertically- 
averaged  current  i;  zero. 

The  initial  u-component  and  v-component  of  velocity  at  100  meters  are  shown  in 
Figure  26  and  Figure  27.  The  u-component  achieves  westward  (dashed  lines)  velocities 
up  to  29  cm  sec-1  in  the  southwestern  portion  of  the  domain.  Weaker  eastward  flow 
up  to  12  cm  sec-1  dominates  the  northern  portion  of  the  domain.  The  v-component 
achieves  southern  (dashed  lines)  speeds  up  to  25  cm  sec-1.  Equatorward  velocities 
dominate  in  the  center  of  the  domain.  Much  weaker  poleward  .  elocitics  up  to  7.5  cm 
sec*1  appear  in  the  northwest  and  southeast  quadrants. 

After  30  hours,  the  u-component  in  the  southwestern  portion  has  intensified  to 
greater  than  40  cm  sec-1  (Figure  28).  Weaker  westward  flow  greater  than  10  cm  sec-1 
is  still  found  in  the  northern  portion  of  the  domain.  The  v-component  has  intensified 
to  greater  than  30  cm  sec-1  in  the  middle  and  southwestern  parts  of  the  domain 
(Figure  29). 
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Cruise  CTZR1,  06-12  July  1988 


T  at  hour  0,  depth  50  m 


Figure  23.  Apparent  Temperature  at  50  m  Depth,  Cruise  CTZRl:  Contour  in¬ 
terval  is  0.2  C. 


The  100  m  pressure  field  (at  time  zero  and  at  time  30  hr),  Figure  30,  shows  what 
mechanism  can  induce  a  jet.  A  high  pressure  center  (dynamic  height  greater  than  15  cm) 
dominates  the  northwestern  portion  of  the  domain,  producing  anticyclonic  flow.  The 
southeastern  portion  of  the  domain  is  dominated  by  a  low  pressure  center  (dynamic 
height  less  than  -15  cm),  producing  cyclonic  flow'.  The  water  between  the  high  and  low 
pressure  centers  is  accelerated  between  the  two  pressure  extremes  which  produces  the 


Figure  24.  Apparent  Temperature  at  100  m  Depth,  Cruise  CTZR1:  Contour  in¬ 
terval  is  0.2  C. 

high  velocities  observed  in  the  jet.  Additionally,  the  strongest  offshore  speeds  corre¬ 
spond  to  the  area  of  strongest  pressure  gradient  between  the  two  pressure  centers  as  seen 
in  the  100  m  velocity  vector  field  at  time  30  hours  in  Figure  31.  The  jet  is  about  50  km 
wide  and  can  be  seen  to  meander  from  the  northeast  to  the  southwest  quadrant  of  the 
domain.  Jessen  and  Ramp  (1988)  reported  speeds  up  to  50  cm  sec-1  across  leg  G  at  the 
194  m  to  210  m  depth  as  .measured  by  ADCP  apparatus. 
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Cruise  CTZR1,  06-12  July  1988 

T  at  hour  0,  depth  200  m 


Figure  25.  Apparent  Temperature  at  200  m  Depth,  Cruise  CTZR1:  Contour  in¬ 
terval  is  0.2  C. 

The  cross-sections  of  the  u-component  and  v-component  velocity  fields  at  time  zero 
are  shown  in  Figure  32  and  Figure  33.  Both  cross-sections  run  northwest  to  southeast 
through  the  area  of  greatest  pressure  gradient  (point  (32.  96)  to  (96,  32)  in  Figure  30). 
Together  the  two  cross-sections  show  the  large-scale  current  patterns  from  the  surface 
to  the  bottom  of  the  ocean.  In  the  u-component  cross-section,  offshore  surface  current 
speeds  are  greater  than  55  cm  sec-1,  decreasing  to  zero  near  300  m.  Below  300  m, 
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denote  offshore  (westward)  flow.  Solid  lines  denote  onshore  (eastward) 
flow.  Contour  interval  is  2.5  cm  sec*1. 


onshore  flow  occurs,  with  speeds  greater  than  15  cm  sec-1  near  700  m.  Flow  gradually 
decreases  to  zero  near  1000  m,  The  v-component  cross-section  shows  equatorward 
surface  current  speeds  greater  than  45  cm  sec*1,  decreasing  to  zero  near  300  m.  Below 
300  m,  poleward  flow  occurs,  with  speeds  greater  than  15  cm  sec*1  near  600  m.  Flow 
again  decreases  to  zero  near  1000  m.  The  large-scale  baroclinic  currents  associated  with 
the  initial  temperature  fields  are  significant  in  the  upper  1000  m  of  the  ocean. 
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Figure  27.  V  Field  at  Time  =  0  h,  Depth  100  m,  Cruise  CTZR1:  Dashed  lines 
denote  equatorward  flow.  Solid  lines  denote  poleward  flow.  Contour 
interval  is  5.0  cm  sec-1. 


B.  VELOCITY  COMPARISON  WITH  ADCP  DATA 

A  qualitative  comparison  between  CTZR1  model  output  velocity  and  ADCP  ac¬ 
quired  velocity  is  now  given.  ADCP  velocities  from  cruise  CTZR1  are  shown  in 
Figure  34.  The  solid  dots  at  regular  intervals  indicate  the  stations  occupied  during  the 
cruise.  Comparison  with  the  model  output  velocity  vector  field  at  time  30  hours, 
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Figure  28.  U  Field  at  Time  =  30  h,  Depth  100  m,  Cruise  CTZR1:  Dashed  lines 
denote  offshore  (westward)  flow.  Solid  lines  denote  onshore  (eastward) 
flow.  Contour  interval  is  5  cm  sec-1. 

Figure  31,  shows  good  agreement  with  the  position  of  the  jet  and  the  magnitude  of  the 
velocities  between  grid  legs  F  and  G. 

Model  domain  cross-section  locations  (solid  lines)  are  shown  in  Figure  35.  A  typi¬ 
cal  cross-section  is  shown  along  grid  leg  C,  Figure  36.  The  corresponding  ADCP 
cross-section  is  shown  in  Figure  37.  In  the  ADCP  cross-sections,  northwest  is  to  the  left 
and  equatorward  flows  are  denoted  by  solid  lines.  In  the  model  output  cross-sections, 
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Figure  29.  V  Field  at  Time  =  30  h,  Depth  100  m,  Cruise  CTZR1:  Dashed  lines 
denote  equatorward  flow.  Solid  lines  denote  poleward  flow.  Contour 
interval  is  5  cm  sec-1. 


northwest  is  to  the  left  and  equatorward  flows  are  denoted  by  dashed  lines.  The  core 
of  the  jet  is  readily  discernable  in  both  cross-sections,  with  the  model  output  showing  a 
60  cm  sec-1  contour.  The  model  output  cross-section  shows  poleward  flow'  (solid  lines) 
in  the  southern  portion  of  the  cross-section,  but  the  ADCP  cross-section  shows  only 
equatorward  flow’  in  the  section  for  comparison.  The  ADCP  cross-section  shows  deeper 
flow's  below'  300  m  (the  irregular  shaded  area  is  the  limit  of  good  data,  not  the  bottom 
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Distonce  (km) 


Figure  30.  Pressure  Field  at  Time  =  0  h,  Depth  100  m,  Cruise  CTZRl:  Contour 
interval  is  2.5  cm.  The  line  segment  joining  points  (32,  96)  and  (96,  32) 
is  the  cross-section  shown  in  Figure  32  and  Figure  33. 

of  the  ocean),  which  indicates  less  shearing  flow  at  300  m  than  the  model  output  cross- 
section,  which  indicates  a  level  of  no  motion  ?.t  300  m.  The  model  output  cross-section 
shows  a  broader  baroclinic  flow  with  a  level  of  no  motion  near  300  m.  In  all  the  cross- 
sections  (except  leg  G),  the  typical  difference  in  surface  velocity  between  theADCP  and 
model  output  velocity  is  5  cm  sec-1  to  10  cm  sec-1,  with  the  ADCP  velocity  usually 
greater  than  the  model  output  velocity. 
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Cruise  CTZR1,  06-12  July  1988 

Vh  at  hour  30,  depth  100  m 


Figure  31.  Velocity  Vector  Field  at  Time  =  30  h,  Depth  100  m,  Cruise 
CTZR1:  Vectors  are  plotted  at  every  other  point  in  the  model  do¬ 
main.  Dots  indicate  positions  of  the  individual  stations  of  the  CTZ 
grid.  Note  the  strongest  flow  between  grid  legs  F  and  G. 


The  cross-sections  across  leg  G  are  of  note  because  of  the  stronger  surface  flows 
shown  in  the  ADCP  velocity  cross-section,  Figure  38,  which  is  about  30  cm  sec-1  higher 
than  the  model  output  cross-section,  shown  in  Figure  39.  The  ADCP  velocity  cross 
section  shows  surface  speeds  greater  than  80  cm  sec-1.  Jessen  and  Ramp  (1988)  reported 
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Cruise  CTZR1,  06-12  July  1988 
u  field  at  hour  0  from  (32,96)  to  (9632) 


NW-SE  Distance  (km) 


Figure  32.  U  Field  Cross-section,  Time  =  0  h,  Depth  100  m,  Cruise 
CTZR1:  Northwest  is  to  the  left.  Westward  flow  is  denoted  by 
dashed  lines.  Eastward  flow  is  denoted  by  solid  lines.  Contour  interval 
is  2.5  cm  sec-1.  Location  of  the  cross-section  is  shown  in  Figure  30. 


surface  speeds  up  to  100  cm  sec-1.  The  model  output  cross-section  only  shows  surface 
speeds  greater  than  55  cm  sec-1. 

It  must  be  remembered  that  the  effect  of  the  objective  analysis  interpolation  between 
20  nm  separated  stations  inevitably  underestimates  geostrophic  flow  with  small  hori¬ 
zontal  scales;  that  is,  tight  gradients  in  the  actual  flow  field  will  be  smoothed  out,  re- 
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Figure  33.  V  Field  Cross-section,  Time  *  0  h,  Depth  100  m,  Cruise 
CTZR1:  Northwest  is  to  the  left.  Equatorward  flow  is  denoted  by 
dashed  lines.  Poleward  flow  is  denoted  by  solid  lines.  Contour  interval 
is  5.0  cm  sec-1.  Location  of  the  cross-section  is  shown  in  Figure  30. 


suiting  in  lower  analyzed  speeds  and  less  horizontal  resolution  in  the  plotted  objectively 
analyzed  flow  field.  This  unavoidable  result  of  the  objective  analysis  is  a  compromise 
between  the  cost  and  resources  available  to  acquire  quasi-synoptic  data  within  a  rea¬ 
sonable  amount  of  time  and  the  need  to  sample  the  flow  field  adequately  to  fully  define 
the  actual  flow'.  Additionally,  the  currents  measured  by  the  ADCP  will  have 
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Figure  34.  ADCP  Velocity  (94-110  m)  for  Cruise  CTZR1:  Regularly  spaced  do*s 
are  stations  occupied  during  the  cruise.  Note  strongest  flows  between 
grid  legs  F  and  G. 
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Figure  35.  Model  Domain  Cross-sections,  Cruise  CTZR1:  Circled  dots  are 
stations  occupied  during  the  cruise.  Solid  lines  show  location  of  model 
cross-sections.  Model  cross-sections  of  legs  C  and  G  will  be  shown  for 
comparison  with  ADCP  data. 


ageostrophic  components  associated  with  them;  e.  g.,  internal  waves,  tides,  inertial  os¬ 
cillations,  and  Ekman  transport. 

All  in  all,  the  model  output  cross-sections  show  a  less  complicated  flow  field  and  a 
shallower  level  of  no  motion  than  the  ADCP  velocity  cross-sections.  To  achieve  a  more 
detailed  and  realistic  deep  flow,  the  barotropic  mode  should  be  incorporated  into  the 
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Figure  36.  Model  Velocity  Cross-section  for  Grid  Leg  C,  Cruise  CTZR1:  Dashed 
lines  indicate  equatorward  flow.  Solid  lines  indicate  poleward  flew. 
Northwest  is  to  the  left.  Contour  interval  is  5  cm  sec-1  . 

model.  A  complete  statistical  verification  of  the  ADCP  acquired  velocity  with  model 
output  velocity  is  considered  to  be  another  study. 

C.  VERTICAL  VELOCITY 

Vertical  velocity  for  CTZR1  at  30  hours  is  shown  in  Figure  40.  The  solid  lines  de¬ 
note  upwelling  areas  and  the  dashed  lines  denote  downwelling  areas.  Strong  mesoscale 
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Figure  37.  ADCP  Velocity  Cross-section  Along  Leg  C  for  Cruise  CTZRl:  Ar¬ 
rows  indicate  boundaries  of  corresponding  model  output  cross-section 
for  comparison.  Solid  lines  indicate  equatorward  flow.  Dashed  lines 
indicate  poleward  flow  Northwest  is  to  the  left.  Contour  interval  is  5 
cm  sec-1. 


patterns  are  clearly  evident.  The  model  calculated  vertical  velocity  by  a  finite  difference 
integration  of  the  continuity  equation  as  given  in  integral  form  by 


*  d*  ,  ,  . 

~ —  dz  —  w(z) 
oz 

The  model  vertical  velocity  value  was  scaled  by  a  factor  of  100  for  plotting  purposes, 
so  the  contour  values  on  the  vertical  velocity  figures  must  be  divided  by  100  to  obtain 
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Figure  38.  ADCP  Velocity  Cross-section  Along  Leg  G  for  Cruise  CTZR1:  Solid 


lines  indicate  equatorvvard  flow.  Northwest  is  to  the  left.  Contour  in¬ 


terval  is  5  cm  sec-1. 


the  values  calculated  by  the  model.  The  maximum  upwelling  speed  was  greater  than  .30 
cm  sec-1’  or  260  m  day1.  The  maximum  downwelling  speed  was  greater  than  .60  cm 
sec-1-  or  520  m  day1.  Washburn  ei  al.,  (1990),  reported  downwelling  rates  of  6  to  10 
m  day'  during  leg  2  (6  to  18  July  1988)  of  the  R/V  THOMAS  WASHINGTON  cruise, 
which  was  coincident  with  the  cruise  dates  for  cruises  CTZR1  and  CTZR2. 
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Figure  39.  Model  Velocity  Cross-sections  for  Leg  G,  Cruise  CTZR1:  Arrows  in¬ 
dicate  boundaries  of  corresponding  ADCP  velocity  cross-section  for 
comparison.  Equatorward  flow  is  denoted  by  dashed  lines.  Poleward 
flow  is  denoted  by  solid  lines.  Northwest  is  to  the  left.  Contour  inter¬ 
val  is  5  cm  sec-1. 


The  areas  of  downwelling  and  upwelling  are  associated  with  the  small  scale  patterns 
of  pressure  (Figure  30).  This  can  be  observed  by  carefully  noting  the  positions  of  the 
areas  of  downwelling  and  upwelling  and  comparing  them  to  the  areas  of  troughing  and 
ridging  in  the  pressure  pattern.  The  areas  of  convergence  are  associated  with  troughs, 
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Figure  40.  Model  Vertical  Velocity  Field  at  Time  30  h,  Depth  100m,  Cruise 
CTZR1:  Downwelling  is  denoted  by  dashed  lines.  Upwelling  is  de¬ 


noted  by  solid  lines.  Contour  interval  is  .10  cm  sec-1. 


and  areas  of  divergence  are  associated  with  ridges.  A  discussion  of  this  relationship  will 
be  pursued  in  the  next  chapter. 

Washburn  et  ai,  (1990)  identified  downwelling  areas  using  bio-optical  data.  The 
points  he  identifies  as  areas  of  downwelling  are  shown  in  Figure  41  and  Figure  42. 

The  major  results  of  this  study  are  contained  in  the  following  three  figures.  The  first 
figure,  Figure  43,  is  of  the  CTZR1  model  vertical  velocity  field  at  30  hours  with  points 
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6-12  JULY,  1988 


PI.  Artno 


Figure  41.  Locations  of  All  Stations  During  Leg  2  of  R/V  THOMAS 
WASHINGTON  Survey  (Triangles):  Filled  triangles  indicate  stations 
with  downwelled  water  masses.  Contours  of  geopotential  anomaly 
(0/500  dbar)  are  from  cruise  CTZR1  (6-12  July  1988).  Offshore  exten¬ 
sion  of  geopotential  anomaly  contours  from  R./V  THOMAS 
WASHINGTON  survey  (from  Washburn  et  al,  1990). 


identified  by  Washburn  et  al,  (1990)  as  downwelling  areas  plotted  onto  the  model  ver¬ 
tical  velocity  field.  Of  the  33  points  plotted,  26  occurred  in  downwelling  areas  (dashed 
lines)  and  seven  points  (filled  in)  occurred  in  upwelling  areas  (solid  lines).  Thus,  an  in¬ 
dependent  estimate  of  downwelling  in  the  CTZ  has  verified  the  model  results,  although 
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Figure  42.  Locations  of  Stations  Containing  Downwelled  Water  Masses:  The  top 
figure  is  from  cruise  CTZR1,  the  middle  figure  is  from  cruise  CTZR2, 
and  the  bottom  figure  is  from  cruise  CTZ2  (from  Washburn  et  al., 
1990). 
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in  a  somewhat  'spotty'  manner.  This  'spottiness'  in  the  verification  charts  for  cruises 
CYZR1,  CTZR2,  and  CTZ2  will  be  addressed  in  the  next  chapter.  The  points  were  from 
the  6  to  12  July  1988  cruises  of  the  R/V  THOMAS  WASHINGTON  and  R/V  POINT 
SUR. 

The  second  figure,  Figure  44,  shows  the  points  identified  from  cruise  CTZR2  as 
downwelling  areas  plotted  on  the  model  vertical  velocity  field  at  time  30  hours.  Of  the 
ten  points  plotted,  eight  points  occurred  in  downwelling  areas,  and  only  two  points  oc¬ 
curred  in  upwelling  areas. 

The  third  figure,  Figure  45,  shows  the  points  identified  from  cruise  CTZ2  as  down¬ 
welling  areas  plotted  on  the  model  vertical  velocity  field  at  30  hours.  Of  the  eight  points 
plotted,  four  occurred  in  downwelling  areas  and  four  occurred  in  upwelling  areas. 

For  all  three  plots,  of  the  51  points  identified  as  downwelling  points,  38  points 
(75%)  occurred  in  downwelling  areas  and  13  occurred  in  upwelling  areas. 

The  magnitude  of  the  mesoscale  vertical  velocity  fields  as  calculated  by  the  model 
were  on  the  order  of  .05  cm  sec-1  for  CTZ2,  .10  cm  sec-1  for  CTZR1,  and  .25  cm  sec-1 
for  CTZR2.  The  unusually  high  order  of  magnitude  of  vertical  velocity  for  cruise 
CTZR2  is  most  probably  a  result  of  the  lack  of  objectively  analyzed  data  below  250  m 
for  this  cruise.  Since  part  of  the  data  processing  involved  fitting  an  exponential  profile 
to  the  last  two  points  of  the  calculated  mean  apparent  temperature  profile,  an  isothermal 
condition  occurs  at  relatively  shallow  depth  (1500  m),  as  seen  in  Figure  46.  The  result¬ 
ing  dynamical  modes  have  large  amplitudes  at  depths  shallower  than  250  m,  and  the 
value  of  the  dynamical  modes  is  essentially  zero  near  1500  m  (Figure  47). 

The  resulting  extrapolation  of  apparent  temperature  to  deep  depths  using  dynamical 
modes  was,  therefore,  forced  by  the  slope  of  the  relatively  short  (  as  compared  to  4000 
m)  250  m  apparent  temperature  profile.  Evidently,  250  m  of  temperature  data  are  not 
sufficient  for  the  accurate  extrapolation  of  temperature  to  depth  using  dynamical  modes. 
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Figure  43.  Verification  of  Model  Vertical  Velocity  Field  at  Time  30  h,  Depth  100m, 
Cruise  CTZR1:  Downwelling  is  denoted  by  dashed  lines.  Upwelling 
is  denoted  by  solid  lines.  Contour  interval  is  0.1  cm  sec*1.  The  points 
identified  as  downwelling  areas  by  Washburn  et  al.  ,  (1990)  and  oc¬ 
curring  in  model  simulated  downwelling  areas  are  circled  with  a  heavy 
line.  Points  identified  as  downwelling  areas  but  occurring  in  model 
simulated  upwelling  areas  are  shaded. 
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Figure  44.  Verification  of  Model  Vertical  Velocity  Field  at  Time  30  h,  Depth  100m, 
Cruise  CTZR2:  Downwelling  is  denoted  by  dashed  lines.  Upwelling 
is  denoted  by  solid  lines.  Contour  interval  is  0.25  cm  sec-1.  The  points 
identified  as  downwelling  areas  by  Washburn  et  al  ,  (1990)  and  oc¬ 
curring  in  model  simulated  downwelling  areas  are  circled  with  a  heavy 
line.  Points  identified  as  downwelling  areas  but  occurring  in  model 
simulated  upwelling  areas  are  shaded. 
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Figure  45.  Verification  of  Model  Vertical  Velocity  Field  at  Time  30  h,  Depth  100m, 
Cruise  CTZ2:  Downwelling  is  denoted  by  dashed  lines.  Upwelling  is 
denoted  by  solid  lines.  Contour  interval  is  0.05  cm  sec-1.  The  points 
identified  as  downwelling  areas  by  Washburn  er  al.  ,  (1990)  and  oc¬ 
curring  in  model  simulated  downwelling  areas  are  circled  with  a  heavy 
line.  Points  identified  as  downwelling  areas  but  occurring  in  model 
simulated  upwelling  areas  are  shaded. 
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Figure  46.  Plot  of  Mean  Apparent  Temperature  Versus  Depth,  Cruise 
CTZR2:  Note  the  isothermal  condition  beginning  near  1500m. 
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V.  DISCUSSION  OF  RESULTS 


Amone  et  al,  (1990)  describes  the  divergence/upwelling  and 
convergence/downwelling  phenomena  as  observed  in  the  Gulf  Stream  system  via  acous¬ 
tic  imaging  of  biological  and  physical  processes  within  Gulf  Stream  meanders.  Evidently 
the  same  divergence/upwelling  and  convergence/downwelling  phenomena  is  present  in 
the  CTZ.  To  gain  some  physical  insight  into  the  relationships  described  in  Amone  et 
al,  (1990),  begin  with  the  two  horizontal  momentum  equations  without  friction: 


and 


du 

dt 


-fv  =  -a 


8P 

6x 


—a 


dP 

dy  ' 


{13} 


{14} 


After  cross-differentiating  ,  subtracting  the  top  equation  from  the  bottom  equation, 
and  doing  appropriate  scaling  analysis,  the  following  simplified  form  of  the  absolute 
vorticity  equation  results: 


{15} 

If  divergence  (V  •  VH  >  0)  is  present,  then  the  absolute  vorticity  tendency  is  negative, 
which  implies  a  tendency  for  anticyclonic  motion.  If  convergence  (V  •  VH  <  0)  is  present, 
ihen  the  absolute  vorticity  tendency  is  positive,  which  implies  a  tendency  for  cyclonic 
motion.  To  aid  in  visualizing  the  preceding  discussion,  Figure  48  is  presented. 

To  conserve  mass,  the  upwelling  centers  are  located  in  areas  of  divergence  situated 
between  the  upstream  trough  and  downstream  ridge.  Similarly,  the  downwelling  centers 
are  located  in  areas  of  convergence  between  the  upstream  ridge  and  downstream  trough. 

Rewriting  •—  (£  +f)  =  -/V  •  V„  as  ~  ({  +J)  =f — -  from  continuity,  then  negative 
absolute  vorticity  tendency  implies  <  0  ,  or  upwelling.  Similarly,  positive  absolute 
vorticity  tendency  implies  ~~  >  0,  or  downwelling. 

The  magnitude  of  upwelling/downwelling  in  the  CTZ  can  be  approximated  by  doing 
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Figure  48.  Convergence  and  Divergence  Patterns  in  a  Typical  Meander:  Stream¬ 
lines  of  a  typical  meander  showing  "the  positions  of  the 
divergence, 'upwelling  centers  and  the  convergence, 'downwelling  centers. 
Vertical  upwelling  centers  are  shown  as  the  fluid  enters  the  crest  in  the 
divergent  region  (dots),  and  downwelling  centers  or  convergence  are 
shown  as  the  fluid  enters  the  trough  (crosses)  (from  Arnone  et  ah, 
1990). 
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a  scale  analysis  of  the/ 


4~  term  in  the  equation  above. 


Now  let 


J  82  H 


{16} 


where 

f  is  the  coriolis  parameter, 

w_H  is  the  vertical  velocity  estimate  at  depth  H,  and 

H  is  the  depth  of  interest  ( 100  m). 

Continuing, 


W-H 

H 


HU 

L 


where  R„  =  — ,  the  Rossby  number. 
J *** 

Based  on  an  initial  estimate  of 


U  ~  10  cm  sec 


the  fluctuating  velocity  averaged  over  the  upper  100  m, 

L  ~  10  Am 


one-fourth  wavelength  of  up'-velling/downwelling  pattern, 

/~  10"4  sec-1 


and 


then 


and 


100  m 


Ro-OA 


{17} 


H 

L 


0.01 
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to  estimate  the  vertical  velocity  as 

w  ~  (.1)(.01)(10)  ~  .Olcw  sec-1 

or  about  10  m  day1,  which  agrees  with  the  upper  estimate  of  Washburn  et  ai,  (1990). 
The  vertical  velocity  magnitudes  from  CTZR1  model  output  are  30  to  60  times  larger 
than  the  above  estimate  from  scale  analysis.  The  cause  of  such  large  differences  between 
estimates  of  vertical  velocity  based  on  observations  and  model  output  is  the  topic  of  a 
future  study. 

The  "spottiness"  of  the  points  identified  as  downwelling  areas  is  now  addressed. 
Washburn  et  ai,  (1990)  estimated  the  rate  of  subduction  occurring  along  the  25.8  and 
26.2  density  surfaces  as  six  to  ten  m  day-'1  and  two  to  four  m  day1,  respectively.  The 
points  shown  by  Washburn  et  ai,  (1990)  in  Figure  41  can  be  considered  as  points 
along  the  subduction  path,  not  isolated  vertical  columns  of  upwelling  or  dow’nwelling. 
At  20  mn  grid  spacing  between  the  alongshore  legs,  exact  points  of  downwelling  and 
upwelling  can  not  be  pin-pointed  by  the  model;  thus,  the  areas  of  downwelling  and  up- 
welling  must  be  considered  as  part  of  larger  mesoscale  features,  not  specific  points  of 
vertical  sinking  or  rising  motion.  A  future  study  should  attempt  to  verify  model  results 
with  observed  downwelling  and  upwelling  areas  at  specific  locations  and  times;  that  is, 
model  runs  should  be  integrated  forward  in  time  until  the  time  of  a  recorded  observation 
(or  set  of  observations)  is  reached,  and  then  the  model  downwelling  and  upwelling  areas 
should  be  compared  to  downwelling  and  upwelling  areas  identified  from  observations. 
In  this  way,  a  more  rigorous  verification  can  be  accomplished.  Acknowledging  the  facts 
that  the  data  points  shown  in  Figure  41  and  Figure  42  (top  figure)  cover  all  six  days 
of  CTZR1,  and  that  the  points  are  not  identified  by  specific  dates  and  times,  it  is  still 
remarkable  that  26  of  33  observed  downwelling  points  occur  on  model  downwelling 
areas.  Similar  statements  applj  to  the  CTZR2  and  CTZ2  data  points. 

This  initial  qualitative  verification  certainly  enlivens  the  possibility  that  a  more  rig¬ 
orous  verification  in  a  future  study  will  show  an  even  higher  percentage  of  observed 
downwelling  or  upwelling  areas  agreeing  with  respective  model  results. 
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VI.  CONCLUSIONS 


The  quality  of  data  acquired  during  cruise  CTZR1  is  sufficiently  high  to  enable  di¬ 
agnosis  of  the  vertical  velocity  field  using  a  numerical  model.  Deeper  CTD  casts  (ideally 
to  the  bottom  of  the  ocean  for  each  station)  would  enable  determination  of  the  first 
three  dynamical  modes  for  use  in  extrapolation  of  temperature  data  to  the  ocean  bot¬ 
tom.  Only  modes  two  and  three  could  be  computed  with  casts  to  500  m.  Attempting 
to  extrapolate  the  apparent  temperature  to  deep  depths  using  data  only  available  to  250 
m  evidently  results  in  unusually  high  vertical  velocity  rates  due  to  an  inaccurate  repre¬ 
sentation  of  the  deep  mean  temperature  profile,  which  is  used  to  calculate  the  dynamical 
modes. 

The  sign  of  the  vertical  velocity  agrees  very  well  with  independent  estimates  from 
Washburn  et  al.,  (1990).  Although  the  model  vertical  velocity  field  and  Washburn's 
data  may  not  be  exactly  in  phase  because  the  data  is  only  quasi-synoptic  in  time,  the 
high  percentage  of  points  (15%)  identified  to  be  in  downwclling  areas  which  occurred  in 
model  output  downwelling  areas  gives  confidence  to  this  verification.,  The  results  are 
certainly  encouraging  and  merit  further  careful  study. 

The  magnitude  of  vertical  velocity  in  the  numerical  model  is  30  to  60  times  larger 
than  the  estimated  vertical  velocity  based  on  observations  from  Washburn  et  al., 
(1990).  This  basic  study  used  a  flat  bottom  with  no  forcing  other  than  the  initial  density 
field  (as  reflected  in  the  apparent  temperature  field  used  to  initialize  the  model).  A  future 
study  with  variable  topography  and  wind  forcing  incorporated  into  the  model  is  neces¬ 
sary  to  determine  their  effect  on  the  magnitude  of  model  output  vertical  velocity. 

The  ADCP  velocity  sections  showed  more  detail  at  deeper  depths  than  the  model 
velocity  sections,  which  showed  broader  flow  at  depth  and  a  level  of  no  motion  near  300 
m.  The  inclusion  of  the  barotropic  mode  into  the  model  should  improve  the  resolution 
of  deep  flows  in  model  velocity  sections.  The  ADCP  and  model  output  velocity  sections 
most  closely  resembled  each  other  in  the  vicinity  of  the  core  of  the  jet.  Differences  be¬ 
tween  the  ADCP  and  model  output  velocity  cross-sections  were  more  pronounced  in 
areas  further  from  the  jet  core. 

The  maximum  surface  velocities  from  model  and  ADCP  cross-sections  agreed  to 
within  10  cm  sec-'  except  for  leg  G,  in  which  the  ADCP  velocity  was  greater  than  the 
model  velocity  by  about  30  cm  sec-1.  Higher  velocities  were  expected  near  leg  G  because 
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of  the  higher  pressure  gradient  in  the  area,  and  the  ADCP  acquired  data  certainly  re¬ 
vealed  this  expectation  to  be  true. 

Future  work  would  include  the  statistical  verification  of  model  velocities  and  ADCP 
velocities,  the  incorporation  of  the  barotropic  mode  into  the  model,  the  use  of  variable 
topography  in  the  model,  and  the  inclusion  of  wind  stress  and  surface  heat  fluxes  in  the 
model.  The  correlation  function  used  in  the  objective  analysis  of  the  CTZ  data  sets 
caused  instabilities  in  the  objectively  analyzed  fields  which  required  a  convective  adjust¬ 
ment  before  initiation  of  the  model  run.  The  determination  of  a  more  optimal  corre¬ 
lation  function  for  the  data  acquired  in  the  CTZ  mapping  cruises  is  necessary  to  fully 
utilize  the  data  sets. 

The  results  of  this  study  are  encouraging  first  steps  in  the  eventual  determination 
of  mass-conserving  velocities  from  numerical  models  initialized  with  quasi-svnoptic  data. 
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